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ABSTRACT 

Using imaging that shows four magnitudes of main sequence stars, we have discovered that the 
Galactic globular cluster NGC 1851 is surrounded by a halo that is visible from the tidal radius of 
700 arcsec (41 pc) to more than 4500 arcsec (>250 pc). This halo is symmetric and falls in density 
as a power law of r — 1 . It contains approximately 0.1% of the dynamical mass of NGC 1851. There 
is no evidence for tidal tails. Current models of globular cluster evolution do not explain this feature, 
although simulations of tidal influences on dwarf spheroidal galaxies qualititivcly mimic these results. 
Given the state of published models it is not possible to decide between creation of this halo from 
isolated cluster evaporation, or from tidal or disk shocking, or from destruction of a dwarf galaxy in 
which this object may have once been embedded. 

Subject headings: Galaxy:Globular Glusters:Individual:NGC 1851, Stars: Hertzsprung-Russell Dia- 
gram 



1. INTRODUCTION 

While it is clear that globular clusters must evapo- 
rate, and that globulars in a tidal field can also lose mass 
through their Lagrangian points, the exact state of dis- 
ruption of many Galactic globulars is less certain. Much 
of the search for tidal structure (and for the fundamental 
parameters of globulars, for that matter) has been done 
with star counts from photographic plates (Grillmair et 
al. 1995, Leon et al. 2000, Illingworth & Illingworth 1976, 
Peterson and King 1975, and references therein). We 
give two examples of problems with the photographic 
datasets: NGC 288, for instance, has been claimed to 
have "extratidal structure" (Leon et al. 2000), while re- 
cent kinematic work finds no evidence for a coherent ve- 
locity gradient in the direction of the tidal tails (Kiss et 
al. 2007). Similarly, Da Costa & Coleman (2008) use 
kinematics to put a limit of 0.7% of the mass of w Cen 
between 1 and 2 tidal radii. 

SDSS, HST, and mosaic cameras on large telescopes 
have allowed new insights into the structure of globu- 
lar clusters. Recently several clusters have been shown 
to have spectacular tidal tails outside the tidal radius: 
Palomar 5 (Odenkirchen et al. 2001, Rockosi et al. 2002, 
Odenkirchen et al. 2003) and NGC 5466 (Belokurov et 
al. 2006, Grillmair & Johnson 2006) are prototypes of this 
class of globulars. In addition, Lauchner et al. (2006) use 
SDSS to claim tidal tails in NGC 5053, though most of 
the stars are confined within the tidal radius. Sohn et 
al. (2003) use CFHT Mosaic imaging to discover "tenta- 
tive detections of tidal halos" around Pal 3 and 4. The 
typical such cluster is "loose" , with a relatively small ra- 
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tio of tidal to core radius. Even those clusters that are 
convincingly being destroyed show the onset of two tidal 
tails at a relatively small distance from the tabulated 
cluster tidal radius. 

It has been known for more than 40 years that the 
main bodies of globular clusters can be described by King 
(1966) models. Limiting (tidal) radii from these King- 
model fits range from 1.6 pc (NGC 6544, galactocentric 
distance 5.3 kpc) to 214 pc (NGC2419, galactocentric 
distance 91.5 kpc), with a median of 36 pc (derived from 
Harris 1996). 90% of all Galactic globulars have radii of 
less than 100 pc (Harris 1996). 

McLaughlin & van der Marel (2005) have shown that 
Wilson (1975) models of the stellar distributions of glob- 
ulars in the Milky Way, the Magellanic Clouds and the 
Fornax dwarf better describe the outer structures of thsc 
clusters than do King models. McLaughlin and van der 
Marel (2005) homogeneously re-fit King, Wilson, and 
power-law models to 153 of these clusters. In general, 
the new King tidal radii are similar to those in the on- 
line Harris catalog, while the Wilson tidal radii are larger. 
Typically, the Wilson tidal radii are a factor of 2.5 larger 
than the King tidal radii (McLaughlin and van der Marel 
2005), while for NGC 1851 the Wilson radius is a factor of 
6.8 larger (McLaughlin and van der Marel's King tidal 
radius for NGC 1851 is 6.7 arcmin as opposed to 11.7 
in the Harris compilation, with a Wilson tidal radius of 
45 armin). We will use the King tidal radii in this pa- 
per simply because they are better known, and since in 
the case of NGC 1851 we believe that the fundamental 
data need improvement (star counts only extend to 10 
arcmin from the center) to determine a realistic limiting 
radius. For the rest of this paper, we will use the phrase 
"tidal radius" to mean the King tidal radius tabulated in 
the online Harris catalog unless stated otherwise. Both 
the King tidal radius and the Wilson tidal radius are 
typically beyond where the star counts end, and can be 
poorly constrained in specific globulars. 

What is the outer structure of a typical globular? Un- 
til recently, NGC 1851, the subject of this paper, would 
have been such a cluster: its absolute magnitude of —8.3 
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is 1 mag brighter than the median Galactic globular, and 
the tidal radius of 41 pc is very close to the the median. 
NGC 1851, however, has recently been shown (Milonc 
et al. 2008), using HST imaging, to be a most remark- 
able object. Above ~mi?606W— 19, there are two dis- 
tinct main-sequences (MS) and lower subgiant branches 
(SGB) visible, making it one of a small number of clusters 
with dual stellar populations. Milonc ct al. (2008) see no 
relative changes in the 2 subgiant populations out to a 
distance of 1.75 arcmin. Zoccali et al. (2009) study the 
radial extent of the two subgiant branches from ground- 
based data, and claim that the fainter sequence dies away 
2.4 arcmin from the center (tidal radius is 11.7 arcmin), 
though Milone et al. (2009) refute the Zoccali work. Di- 
nescu et al. (1997) derive the absolute proper motion for 
NGC 1851, which when coupled with the large (320.5 
km/s) heliocentric radial velocity and current distances 
(R Q =12.1 kpc, R ffC =16.7 kpc, X= -4.3 kpc, Y= -8.9 
kpc, Z= —6.9 kpc), gives an orbit with period ~0.4 Gyr, 
apogalacticon distance ^30 kpc, and perigalacticon dis- 
tance ~ 5 kpc. This orbit passes through the disk of the 
Milky Way five times per Gyr at distances ranging from 
5 to 30 kpc (orbital calculations graciously performed 
by Piatek (2009)). According to the theoretical papers 
cited in the discussion, disk shocking and tidal shock- 
ing should be important for NGC 1851, and a tidal tail 
should be seen. Leon et al (2000, their figure 10) claim 
to see symmetric structure to about 20 arcmin radius 
(about 1.7x the tidal radius) and then "extensions which 
are likely tracing the orbital path..." Leon et al (2000) 
azimuthally average background-count subtracted data 
and fit a power law to radially distant points, finding a 
power law slope for NGC 1851 of —0.98, similar to what 
we find in Section 2. Our new observations contradict 
these claims: we see no tidal tails and do not see discrete 
structures in common with the photographic counts, but 
show smooth extended structure. 

2. OBSERVATIONS AND DATA REDUCTION 

We discovered that NGC 1851 has extended structure 
while observing a control field, centered 3.5 tidal radii 
from NGC 1851, for a project to probe the outer struc- 
ture of the LMC. The R versus C-R (C is Washington C 
(Canterna 1976), centered at ^4000 A) color-magnitude 
diagram (CMD) of this field C18 (1=245, b= -35) is com- 
pared to C20 (1=245, b= -25), 10 degrees away, in Fig- 
ure [TJ A main sequence is obvious in C18 from R~18.5 
to R~22. C20, at lower latitude, has more stars overall, 
but does not show this additional feature. This feature is 
most likely a main sequence of an object at a single dis- 
tance, and corresponds exactly to the main-sequence lo- 
cus in published NGC 1851 main body photometry (e.g., 
Figure 10 in Walker 1992), after a reasonable color trans- 
formation. 

We therefore began to map out this main sequence by 
acquiring nine fields outlining the cluster from ~1.0-~6.5 
tidal radii. Exposures were 300 seconds in V and 600 
seconds in B, using the CTIO 4m Blanco telescope with 
the Mosaic II CCD array. Each field is 36x36 arcmin on 
a side. Table Q] gives a diary of observations of the BV 
fields. 

Data were reduced using IRAF's 7 MSCRED package 

7 IRAF is distributed by the National Optical Astronomy Ob- 



and a set of homegrown IDL routines (Saha et al. (2010). 
Photometry is then done with a privately modified (by 
Saha) version of DoPHOT (Schechter et al. 1993) and 
with a set of post-DoPHOT programs written by one of 
us (Saha). Two aspects of the reduction are relatively un- 
common: first, we create single resampled images from 
the 8-CCD mosaics as input to DoPHOT. Our exten- 
sive testing shows that the photometry is unaffected by 
our resampling of the images (Saha et al. 2010). Second, 
DoPHOT was run with a non- variable PSF. However, the 
real PSF does vary with position on the field-of-view, so 
the PSF fitted magnitudes have consequent position de- 
pendent trends of several percent. These were corrected 
by obtaining aperture magnitudes for the brighter high 
S/N stars (which do not vary with position), and compar- 
ing the PSF fitted magnitudes to the aperture measured 
values, and thus deducing the position dependence of the 
fitted magnitudes. The derived corrections as a function 
of position on the field-of-view were then applied for all 
stars. This procedure is more robust for these data than 
using a variable PSF. These extra steps also allow for cor- 
rections due to tilt and piston and defocus as a function 
of position on each original CCD. There is demonstrably 
little chip-to-chip photometry shift remaining after using 
this technique. The small chip-to-chip and field-to-field 
photometry errors from these reduction procedures keep 
the photometry from wandering in zeropoint and scale, 
and permit the large-scale mapping described here. 

The nine fields are combined by using stars in common 
between frames to adjust the instrumental magnitudes 
to a common scale, the scale of our innermost field. The 
photometric zeropoint of the innermost field comes from 
comparison of published NGC 1851 photometry (Walker 
1992). In other words, photometry within ~1 tidal ra- 
dius was forced to that of Walker (1992), but the main 
sequence outside that radius is bootstrapped to the mag- 
nitudes so obtained in the innermost field (field 2 shifted 
by using overlap region of fields 1 and 2, field 3 shifted 
by using overlap of fields 1 and 2 against it, etc.) Thus 
any excursions of main sequence structure over the ex- 
tended field are preserved, if they exist. The upper plot 
of Figure [2] shows the CMD of the entire area imaged, 
including the small region inside 1 tidal radius. We see a 
clear main sequence and lower subgiant branch, the wall 
of thick disk and halo turnoff stars, the wall of K and 
M stars, and the blob of faint galaxies. The lower plot 
shows only stars more than 2000 arcsec from the cluster 
center (more than 2.85 tidal radii). In both CMDs there 
are two turnoffs and lower subgiant branches (Figure [3]), 
extending the Milone et al. (2008) HST result to extraor- 
dinarily large radii. This result contradicts the Zoccali 
et al. (2009) claim that the fainter SBC disappears 2.4 
arcmin from the cluster center, and is consistent with the 
recent Milone et al. (2009) result which follows the two 
subgiant branches out to at least 8 arcmin. 

We define a region on the CMD that contains the main 
sequence and lower subgiant branch. This region is a 
band that follows the ridge line of the MS and SGB and 
is 0.15 mag wide in (B-V) (see Figure 01. Background 
estimates come from regions centered 0.15 mag redward 

servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. — R versus C-R cmds for control fields C20 and C18 (near NGC 1851) from another project. Note the wall of F stars and the 
"extra" main sequence in C18. 



TABLE 1 

Log of BV observations of NGC 1851 



Field RA J2000.0 Dec J2000.0 M.J.D. of V data 



NGC 1851 inner 
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15 


06.8 


-40 


26 


47 


54825.20172858 


NGC 1851 E 


05 


16 


57.1 


-40 


02 


51 


54715.40055975 


NGC 1851 W 


05 


11 


16.0 


-40 


03 


31 


54823.28117295 


NGC 1851 N 


05 


14 


07.1 


-39 


29 


57 


54822.33423601 


NGC 1851 S 


05 


14 


07.2 


-40 


35 


50 


54715.38609743 


NGC 1851 NE 
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16 


59.9 
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15.0 
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-40 


35 
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54823.29429358 
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35 


51 


54823.30709821 
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Fig. 2. — V versus B-V cmds for two regions surrounding NGC 1851. Top diagram is of all acquired data. Bottom figure is of a region 
more than 2.85 tidal radii from NGC 1851. In both cases the main sequence is identical to that of the globular, and the double subgiant 
branches are visible. 
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Fig. 3. — Histograms of all stars with 19.0 < V < 19.4, in B-V bins of 0.02, to show peaks in the CMD. The features at B-V=0.55 
and 0.63 are the two stellar populations in NGC 1851. The peak at about B— V=1.4 is made up of the foreground K and M stars. Top: 
Histogram of all observed stars; Middle: Histogram of all stars closer than 2000 arcsec (117 pc); Bottom: Histogram of all stars more than 
2400 arcsec (140 pc) from the center to about 4000 arcsec (230 pc). The two peaks of interest are marked with tickmarks in the bottom 
plot. 
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and blueward of this box, and by scaling the C20 control 
field in Figure [TJ In Figure [5] we plot the positions of 
stars in our MS and SGB regions readily distinguished 
in Figure [2j We also plot the distribution of stars in 
the two background regions of the CMD. Aside from a 
density peak of non-cluster stars at large distances to the 
East, there are no obvious trends. Figure [5] shows the 
angular distribution of stars outside the tidal radius, in 
15° bins, of stars in the main-sequence region of the CMD 
(dashed), in the two background regions of the CMD 
(dotted), and in the total number of NGC stars (solid), 
which is the counts in the dashed histogram minus half of 
the dotted. Counting errors for the solid curve are ±7 on 
average, which encompasses most of the bins. There is 
no evidence, therefore, for tidal tails or angular structure 
on the distribution of NGC 1851 stars. 

Counting stars on the main sequence, and counting 
stars to the red and blue of the main sequence to use 
as background estimations, gives the density plot of Fig- 
ure [7] We see steeply declining counts from about 500 
arcsec to about 1000 arcsec (1.3 x the cataloged tidal ra- 
dius), with a very shallow decline out to the limits of our 
imaging, more than 4000 arcsec. Fitting a straight line 
to this halo, beyond 1200 arcsec, gives log(counts) oc 
log(r) 24 , with 1-a errors on the power law of —0.58 
and —1.9. Figure [8] shows the density plot in linear 
units to underscore the large extent of the NGC 185 halo. 
While our discovery of stars out to 67 arcmin may not 
seem extreme if one believes the Wilson tidal radius of 45 
arcmin for NGC 1851, we note that the King model fit 
using the code of McLaughlin and van der Marel (2005) 
gives no stars beyond 5.8 arcmin, and the Wilson model 
gives no stars beyond 43 arcmin. We can still see the 
main sequence by eye (Figure [2|) at radii extending be- 
yond the (very uncertain)Wilson tidal radius. 

3. DISCUSSION 

Unlike the case of the Pal 5 tidal tails, this set of ex- 
tratidal stars does not make up a large fraction of the 
current mass of the cluster. For a region exterior to 
13 arcmin from the center, we count 1647 total stars 
in the main sequence region of Figure 2J Using off-main- 
sequence regions defined in Figure [4] and in a scaled ver- 
sion of control field C20, we measure that 699 of those 
stars are field stars. The net number of stars outside 13 
arcmin is therefore 948±40.5, or 58% of the total. De- 
spite the extremely low surface brightness of ^32 mag 
arcsec -2 , the main sequence and the SGB remain vivid. 

If the typical mass on this main sequence is 0.5 M , 
this halo of 500 M Q measures 0.1% of the dynamical mass 
of NGC 1851. NGC 1851 is therefore embedded in a very 
low-mass halo 500pc across, or "dwarf spheroidal size" . 
This halo also contains two populations, as seen in the 
two main sequences, of unknown origin. 

The state of modelling of stars outside the King tidal 
radius and not in tidal tails, is unsatisfying. Observa- 
tionally, as pointed out in Odenkirchen et al. (2009), the 
expected signatures of tidal mass loss are not as clear 
as hoped. Gnedin and Ostriker (1997) model the de- 
struction rates of globulars, but do not provide detailed 
spatial stellar distributions. A number of papers (e.g., 
Spitzer 1940; Gnedin & Ostriker 1997) show that mass 
loss, because of internal 2-body interactions that force 
lower- mass stars slowly to the periphery of the cluster, 



changes the mass- function from center to edge of the clus- 
ter. Combes et al. (1999) find that unbound stars (not 
in the tidal tails) change slope at the tidal radius, to a 
surface density slope of —3. Combes et al. show a time 
sequence of their model m2 (their table 1 and figure 14). 
m2 has a tidal radius of 60 pc. Even after significant 
mass loss the inner symmetrical part stays at about 60 
pc radius, with the tidal tails emanating from the tidal 
radius area. Fukushige & Heggie (2000) use simulations 
to show that stars inside the tidal radius with energies 
above the escape energy can leave the cluster on a similar 
timescale to stars leaving due to 2-body interactions, ap- 
proximately a Hubble time. This simulation might come 
closest to mimicking a cluster like NGC 1851, except for 
the complications due to the orbital eccentricity of NGC 
1851. Mass loss in a slightly different context, that of 
dwarf spheroidal galaxies, has been studied by Johnston 
et al. (1999) and by Peharrubia et al. (2008). Johnston 
et al (1999), in a study of the mass loss rates from dwarf 
spheroidal galaxies, point out the existance of a change 
of slope of the stellar density (star counts) at the radius 
where now- unbound stars become important. They de- 
rive a much shallower slope for the region outside the 
tidal radius, —1 or so, in contrast to that of Combes et 
al. 

Typical observational cases might be those of NGC 
2419 and M92. Ripepi et al. (2007), see main sequence 
stars in NGC 2419 out to more than 1.2 tidal radii. Lee 
et al. (2003) show the radial profile for M92, with power 
law fits of slope ~ — 1 beyond the tidal radius, but with 
symmetric structure out to perhaps 1.5 x the tidal ra- 
dius, little evidence for tidal tails, and possible clumpy 
structure outside the symmetric structure. This is quali- 
tatively consistent with Testa et al. (2000) whose analysis 
of POSS data shows symmetric structure out to perhaps 
1.5-2 x the tidal radius, plus more distant irregular struc- 
ture. Testa et al. counted stars to 16th magnitude in 
photographic F, while Lee et al. counted stars to V=20.0 
and to 23.5. Some of the radially distant structures are in 
common to both studies, at least qualitatively. However, 
Drukier et al. (2007) show that beyond the M92 tidal ra- 
dius, almost no stars have the M92 velocity (their Figure 

In addition, good evidence for extratidal material in 
some globulars has been overlooked. Lee et al. (2004, 
their Figure 5) show the CMD of a control field at more 
than 2.5 tidal radii from NGC 7492 that clearly shows 
the cluster main sequence. 

Given the state of the modelling described above, it is 
impossible to distinguish between creation of such a halo 
from 2-body internal interactions, disk shock heating, 
tidal heating, destruction of a primordial dwarf galaxy in 
which NGC 1851 was embedded, and initial conditions. 
Nonetheless, we briefly examine some of these. 

Mass loss caused by Galactic tides is the obvious ex- 
planation for the long tidal tails of objects such as 
Pal 5 (Odenkirchen et al. 2001; Rockosi et al. 2002; 
Odenkirchen et al. 2003; Koch et al. 2004). The Pal 
5 tidal tails emerge from the cluster at about the tidal 
radius of 16 arcmin, or ^100 pc (Odenkirchen et al. 2001, 
their Figure 2). We see no sign of tidal tails in the NGC 
1851 halo out to the limit of our survey, a radius of 250 
pc. The orbit of Pal 5 (Scholz et al. 1998), with R m i n 
between 5 and 9 kpc, R m ax between 17 and 18 kpc, and 
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Fig. 4. — CMD of NGC 1851 halo with our main-sequence and subgiant regions, and the two off-sequence background regions overlaid. 
Each region is 0.15 mag wide in B-V. There is a very small overlap at V=18.5. 
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Fig. 5. — Left: positions of all stars that meet our criteria for being on the observed main sequence and subgiant branches. Regions 
interior to 1 tidal radius (700 arcsec) were not observed except for one small region to the south. The extreme plotted points are more 
than 6 tidal radii from the center. Right: the same plot for the distribution of stars in the two background regions. An excess of field stars 
can be seen to the East. 
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Fig. 6. — The angular distribution of stars as a proof of the lack of tidal tails. Dashed histogram: all stars outside of ~1 tidal radius, 
in the main sequence region of Figure [4] Dotted histogram: all stars in the two off-cluster positions (twice the area in the CMD) shown 
in Figure [4] Solid histogram: the difference "dashed" minus 0.5x"dotted" (there are two off-main-sequence regions). This is the total 
number of NGC 1851 stars in the main sequence region. The average is about 25±7 for the solid histogram, showing the lack of statistically 
significant features in the azimuthal distribution. 
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Fig. 7. — Density plot of stars surrounding NGC 1851. Note the Y scale of log of stars per 64 arcmin 2 .Open circles are regions in radius 
with very incomplete coverage. Filled circles between 1200 and 4000 arcsec were used for power-law fits. Solid line is the best fit (x 2 = 0.9) 
of r -1 - 24 ; dashed lines are the l-<r limits of r~ 19 and r~ - 58 . 



orbital period 0.4 Gyr, is qualitatively similar to that de- 
rived for NGC 1851 above. Since Pal 5 is much fluffier 
than NGC 1851, we would indeed expect that Pal 5 would 
be closer to dissolution than NGC 1851, as it is. Leon 
et al. (2000) claim evidence for tidal tails tracing NGC 
1851's orbit, evident at distances of 15-40 arcmin (50-140 
pc). Our CMDs contradict these claims. 

Gnedin and Ostrikcr (1997) show that the evaporation 
time for NGC 1851 is very long, about 10 Hubble times. 
Gravitational shocking is also of order the same value. If 
the halo of NGC 1851 is simply an evaporated or shocked 
halo, there is qualitative consistency with Pal 5: Pal 5 
(from their table 3) is quickly destroyed, while NGC 1851 
is very slowly destroyed. The models do not provide the 
information necessary to compare our empirical halo with 
a model halo. 

The studies of destruction of dwarf spheroidals (John- 
ston et al. 1999; Penarrubia et al. 2008, their Figure 2) of- 
fer a possible explanation for the structure of NGC 1851. 
The Penarrubia et al. paper models a dwarf spheroidal 
as an NFW dark-matter profile with a much more con- 
centrated baryonic component. NGC 1851 would be 
an even-more concentrated nucleus in such a model. 
Penarrubia et al. (2008, their figures 2 and 3) show the 
effect of stripping away most of the mass of the dSph. 
What remains is an object that can be fit to a King 
model, but with a shallow halo outside the place where 
the crossing time equals half the orbital time. One prob- 
ably cannot extend these models to the remnants of a 
nucleated dwarf, so we must appeal to n-body experts 
to make the proper models for NGC 1851. The most 
striking feature of NGC 1851 when compared to ui Cen, 



M54, and NGC 2419, clusters that have at one time been 
hoped to be remnant nuclei of nucleated dwarfs, is the 
small core radius, 0.2 pc, of NGC 1851. NGC 1851 (Yong 
et al. 2009, Ventura et al. 2009 ), similarly to lo Cen, 
does have a large spread in the elemental abundances 
C+N+O, making it different from most globulars. 

Of the 24 Galactic globulars more luminous than NGC 
1851, the four clusters with the largest current galacto- 
centric distances all have larger Wilson tidal radii than 
NGC 1851, and three have interesting associations with 
stellar streams: NGC 1851 is on an elliptical orbit reach- 
ing in to a galactocentric distance of 5 kpc and out to 
30 kpc (Section 1). NGC 2419 might be associated with 
the Virgo Stellar Stream, which is on an extremely el- 
liptical orbit (Casctti-Dinescu et al. 2009). NGC 5824 
may be associated with a stellar stream (Newberg et 
al. 2009). NGC 6715 is M54, the core of the Sgr dSph. 
The other three luminous clusters with larger Wilson 
tidal radii than NGC 1851 are interior to the current 
position of NGC 1851, with NGC 6139 the most extreme 
(R GC = 1-7 kpc and Rtidal{Wilson) = 676 pc.) The orbit 
of 47 Tuc is rather circular in the XY plane (R apo = 7.3 
kpc, Rp er i = 5.2kpc), while M3's orbit has an cllipticity 
of 0.42 (R apo = 13.4 kpc, R per i = 5.5kpc)(Dinescu et 
al. 1999). A parameter that relates luminous globulars 
to large physical size might therefore be some combina- 
tion of short-period orbit, eccentric orbit, or association 
with a stellar stream. 

While the correct explanation for the halo of NGC 1851 
is presently elusive, the observational imaging result is 
clear and strong. Modern imaging and reduction tech- 
niques can therefore tell if NGC 1851 is an exceptional 



10 



c 

*E 
i 

u 



-i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i — i — r- 



CT 

if) 

CD 

1_ 
CD 
CL 

Cfl 

o 

H — ' 

w 
CO 



1 - 



J. 



s • 

•A .V* • *•*.: •* . : 

: . • • • • I 



0~) 
00 

O 
O 



- 



cx> _ 1 i- 
o 1 i 







_l I I I I I I I I I I I I I I I I I I I I I I I fi_l I I I I I I I I I I I I I I I L 

1000 2000 3000 4000 5000 

Dist from cluster (orc-sec) 



Fig. 8. — Density plot of stars surrounding NGC 1851. The x axis is now linear in radius. 



cluster and if the fundamental data for the 150 Milky 
Way globular clusters need to be much improved and 
extended with wide-field CCD photometry. 

We thank the CTIO mountain staff for their help in 
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